Evaluating the utility of a high throughput thiol-containing fluorescent probe to screen for reactivity: A case study with the Tox21 library.

High-throughput screening (HTS) assays for bioactivity in the Tox21 program aim to evaluate an array of different biological targets and pathways, but a significant barrier to interpretation of these data is the lack of high-throughput screening (HTS) assays intended to identify non-specific reactive chemicals. This is an important aspect for prioritising chemicals to test in specific assays, identifying promiscuous chemicals based on their reactivity, as well as addressing hazards such as skin sensitisation which are not necessarily initiated by a receptor-mediated effect but act through a non-specific mechanism. Herein, a fluorescence-based HTS assay that allows the identification of thiol-reactive compounds was used to screen 7,872 unique chemicals in the Tox21 10K chemical library. Active chemicals were compared with profiling outcomes using structural alerts encoding electrophilic information. Random Forest classification models based on chemical fingerprints were developed to predict assay outcomes and evaluated through 10-fold stratified cross validation (CV). The mean CV Balanced Accuracy of the validation set was 0.648. The model developed shows promise as a tool to screen untested chemicals for their potential electrophilic reactivity based solely on chemical structural features.

Verification of In Vivo Estrogenic Activity for Four Per- and Polyfluoroalkyl Substances (PFAS) Identified as Estrogen Receptor Agonists via New Approach Methodologies.

Given concerns about potential toxicological hazards of the thousands of data-poor per- and polyfluorinated alkyl substances (PFAS) currently in commerce and detected in the environment, tiered testing strategies that employ high-throughput in vitro screening as an initial testing tier have been implemented. The present study evaluated the effectiveness of previous in vitro screening for identifying PFAS capable, or incapable, of inducing estrogenic responses in fish exposed in vivo. Fathead minnows (Pimephales promelas) were exposed for 96 h to five PFAS (perfluorooctanoic acid [PFOA]; 1H,1H,8H,8H-perfluorooctane-1,8-diol [FC8-diol]; 1H,1H,10H,10H-perfluorodecane-1,10-diol [FC10-diol]; 1H,1H,8H,8H-perfluoro-3,6-dioxaoctane-1,8-diol [FC8-DOD]; and perfluoro-2-methyl-3-oxahexanoic acid [HFPO-DA]) that showed varying levels of in vitro estrogenic potency. In agreement with in vitro screening results, exposure to FC8-diol, FC10-diol, and FC8-DOD caused concentration-dependent increases in the expression of transcript coding for vitellogenin and estrogen receptor alpha and reduced expression of insulin-like growth factor and apolipoprotein eb. Once differences in bioconcentration were accounted for, the rank order of potency in vivo matched that determined in vitro. These results provide a screening level benchmark for worst-case estimates of potential estrogenic hazards of PFAS and a basis for identifying structurally similar PFAS to scrutinize for putative estrogenic activity.

Evaluation of 147 perfluoroalkyl substances for immunotoxic and other (patho)physiological activities through phenotypic screening of human primary cells.

A structurally diverse set of 147 per- and polyfluoroalkyl substances (PFAS) was screened in a panel of 12 human primary cell systems by measuring 148 biomarkers relevant to (patho)physiological pathways to inform hypotheses about potential mechanistic effects of data-poor PFAS in human model systems. This analysis focused on immunosuppressive activity, which was previously reported as an in vivo effect of perfluorooctanoic acid (PFOA) and perfluorooctanesulfonic acid (PFOS), by comparing PFAS responses to four pharmacological immunosuppressants. The PFOS response profile had little correlation with reference immunosuppressants, suggesting in vivo activity does not occur by similar mechanisms. The PFOA response profile did share features with the profile of dexamethasone, although some distinct features were lacking. Other PFAS, including 2,2,3,3-tetrafluoropropyl acrylate, demonstrated more similarity to the reference immunosuppressants but with additional activities not found in the reference immunosuppressive drugs. Correlation of PFAS profiles with a database of environmental chemical responses and pharmacological probes identified potential mechanisms of bioactivity for some PFAS, including responses similar to ubiquitin ligase inhibitors, deubiquitylating enzyme (DUB) inhibitors, and thioredoxin reductase inhibitors. Approximately 21% of the 147 PFAS with confirmed sample quality were bioactive at nominal testing concentrations in the 1-60 micromolar range in these human primary cell systems. These data provide new hypotheses for mechanisms of action for a subset of PFAS and may further aid in development of a PFAS categorization strategy useful in safety assessment.

Comprehensive interpretation of in vitro micronucleus test results for 292 chemicals: from hazard identification to risk assessment application.

Risk assessments are increasingly reliant on information from in vitro assays. The in vitro micronucleus test (MNvit) is a genotoxicity test that detects chromosomal abnormalities, including chromosome breakage (clastogenicity) and/or whole chromosome loss (aneugenicity). In this study, MNvit datasets for 292 chemicals, generated by the US EPA's ToxCast program, were evaluated using a decision tree-based pipeline for hazard identification. Chemicals were tested with 19 concentrations (n = 1) up to 200 µM, in the presence and absence of Aroclor 1254-induced rat liver S9. To identify clastogenic chemicals, %MN values at each concentration were compared to a distribution of batch-specific solvent controls; this was followed by cytotoxicity assessment and benchmark concentration (BMC) analyses. The approach classified 157 substances as positives, 25 as negatives, and 110 as inconclusive. Using the approach described in Bryce et al. (Environ Mol Mutagen 52:280-286, 2011), we identified 15 (5%) aneugens. IVIVE (in vitro to in vivo extrapolation) was employed to convert BMCs into administered equivalent doses (AEDs). Where possible, AEDs were compared to points of departure (PODs) for traditional genotoxicity endpoints; AEDs were generally lower than PODs based on in vivo endpoints. To facilitate interpretation of in vitro MN assay concentration-response data for risk assessment, exposure estimates were utilized to calculate bioactivity exposure ratio (BER) values. BERs for 50 clastogens and two aneugens had AEDs that approached exposure estimates (i.e., BER < 100); these chemicals might be considered priorities for additional testing. This work provides a framework for the use of high-throughput in vitro genotoxicity testing for priority setting and chemical risk assessment.

A gene expression biomarker for predictive toxicology to identify chemical modulators of NF-κB.

The nuclear factor-kappa B (NF-κB) is a transcription factor with important roles in inflammation, immune response, and oncogenesis. Dysregulation of NF-κB signaling is associated with inflammation and certain cancers. We developed a gene expression biomarker predictive of NF-κB modulation and used the biomarker to screen a large compendia of gene expression data. The biomarker consists of 108 genes responsive to tumor necrosis factor α in the absence but not the presence of IκB, an inhibitor of NF-κB. Using a set of 450 profiles from cells treated with immunomodulatory factors with known NF-κB activity, the balanced accuracy for prediction of NF-κB activation was > 90%. The biomarker was used to screen a microarray compendium consisting of 12,061 microarray comparisons from human cells exposed to 2,672 individual chemicals to identify chemicals that could cause toxic effects through NF-κB. There were 215 and 49 chemicals that were identified as putative or known NF-κB activators or suppressors, respectively. NF-κB activators were also identified using two high-throughput screening assays; 165 out of the ~3,800 chemicals (ToxCast assay) and 55 out of ~7,500 unique compounds (Tox21 assay) were identified as potential activators. A set of 32 chemicals not previously associated with NF-κB activation and which partially overlapped between the different screens were selected for validation in wild-type and NFKB1-null HeLa cells. Using RT-qPCR and targeted RNA-Seq, 31 of the 32 chemicals were confirmed to be NF-κB activators. These results comprehensively identify a set of chemicals that could cause toxic effects through NF-κB.

Comprehensive assessment of NR ligand polypharmacology by a multiplex reporter NR assay.

Nuclear receptors (NR) are ligand-modulated transcription factors that regulate multiple cell functions and thus represent excellent drug targets. However, due to a considerable NR structural homology, NR ligands often interact with multiple receptors. Here, we describe a multiplex reporter assay (the FACTORIAL NR) that enables parallel assessment of NR ligand activity across all 48 human NRs. The assay comprises one-hybrid GAL4-NR reporter modules transiently transfected into test cells. To evaluate the reporter activity, we assessed their RNA transcripts. We used a homogeneous RNA detection approach that afforded equal detection efficacy and permitted the multiplex detection in a single-well format. For validation, we examined a panel of selective NR ligands and polypharmacological agonists and antagonists of the progestin, estrogen, PPAR, ERR, and ROR receptors. The assay produced highly reproducible NR activity profiles (r > 0.96) permitting quantitative assessment of individual NR responses. The inferred EC50 values agreed with the published data. The assay showed excellent quality ( = 0.73) and low variability ( = 7.2%). Furthermore, the assay permitted distinguishing direct and non-direct NR responses to ligands. Therefore, the FACTORIAL NR enables comprehensive evaluation of NR ligand polypharmacology.

A Roadmap to the Structure-Related Metabolism Pathways of Per- and Polyfluoroalkyl Substances in the Early Life Stages of Zebrafish (Danio rerio).

BACKGROUND: Thousands of per- and polyfluoroalkyl substances (PFAS) with diverse structures have been detected in the ambient environment. Apart from a few well-studied PFAS, the structure-related toxicokinetics of a broader set of PFAS remain unclear. OBJECTIVES: To understand the toxicokinetics of PFAS, we attempted to characterize the metabolism pathways of 74 structurally diverse PFAS samples from the U.S. Environmental Protection Agency's PFAS screening library. METHODS: Using the early life stages of zebrafish (Danio rerio) as a model, we determined the bioconcentration factors and phenotypic toxicities of 74 PFAS. Then, we applied high-resolution mass spectrometry-based nontargeted analysis to identify metabolites of PFAS in zebrafish larvae after 5 d of exposure by incorporating retention time and mass spectra. In vitro enzymatic activity experiments with human recombinant liver carboxylesterase (hCES1) were employed to validate the structure-related hydrolysis of 11 selected PFAS. RESULTS: Our findings identified five structural categories of PFAS prone to metabolism. The metabolism pathways of PFAS were highly related to their structures as exemplified by fluorotelomer alcohols that the predominance of ß-oxidation or taurine conjugation pathways were primarily determined by the number of hydrocarbons. Hydrolysis was identified as a major metabolism pathway for diverse PFAS, and perfluoroalkyl carboxamides showed the highest in vivo hydrolysis rates, followed by carboxyesters and sulfonamides. The hydrolysis of PFAS was verified with recombinant hCES1, with strong substrate preferences toward perfluoroalkyl carboxamides. CONCLUSIONS: We suggest that the roadmap of the structure-related metabolism pathways of PFAS established in this study would provide a starting point to inform the potential health risks of other PFAS. https://doi.org/10.1289/EHP7169.

Quantitative Chemical Proteomics Reveals Interspecies Variations on Binding Schemes of L-FABP with Perfluorooctanesulfonate.

Evaluating interspecies toxicity variation is a long-standing challenge for chemical hazard assessment. This study developed a quantitative interspecies thermal shift assay (QITSA) for in situ, quantitative, and modest-throughput investigation of chemical-protein interactions in cell and tissue samples across species. By using liver fatty acid binding protein (L-FABP) as a case study, the QITSA method was benchmarked with six per- and polyfluoroalkyl substances, and thermal shifts (ΔTm) were inversely related to their dissociation constants (R2 = 0.98). The QITSA can also distinguish binding modes of chemicals exemplified by palmitic acid. The QITSA was applied to determine the interactions between perfluorooctanesulfonate (PFOS) and L-FABP in liver cells or tissues from humans, mice, rats, and zebrafish. The largest thermal stability enhancement by PFOS was observed for human L-FABP followed by the mouse, rat, and zebrafish. While endogenous ligands were revealed to partially contribute to the large interspecies variation, recombinant proteins were employed to confirm the high binding affinity of PFOS to human L-FABP, compared to the rat and mouse. This study implemented an experimental strategy to characterize chemical-protein interactions across species, and future application of QITSA to other chemical contaminants is of great interest.

Bioactivity profiling of per- and polyfluoroalkyl substances (PFAS) identifies potential toxicity pathways related to molecular structure.

Per- and polyfluoroalkyl substances (PFAS) are a broad class of hundreds of fluorinated chemicals with environmental health concerns due to their widespread presence and persistence in the environment. Several of these chemicals have been comprehensively studied for experimental toxicity, environmental fate and exposure, and human epidemiology; however, most chemicals have limited or no data available. To inform methods for prioritizing these data-poor chemicals for detailed toxicity studies, we evaluated 142 PFAS using an in vitro screening platform consisting of two multiplexed transactivation assays encompassing 81 diverse transcription factor activities and tested in concentration-response format ranging from 137 nM to 300 µM. Results showed activity for various nuclear receptors, including three known PFAS targets--specifically estrogen receptor alpha and peroxisome proliferator receptors alpha and gamma. We also report activity against the retinoid X receptor beta, the key heterodimeric partner of type II, non-steroidal nuclear receptors. Additional activities were found against the pregnane X receptor, nuclear receptor related-1 protein, and nuclear factor erythroid 2-related factor 2, a sensor of oxidative stress. Using orthogonal assay approaches, we confirmed activity of representative PFAS against several of these targets. Finally, we identified key PFAS structural features associated with nuclear receptor activity that can inform future predictive models for use in prioritizing chemicals for risk assessment and in the design of new structures devoid of biological activity.

Characterisation and validation of an in vitro transactivation assay based on the 22Rv1/MMTV_GR-KO cell line to detect human androgen receptor agonists and antagonists.

We describe the characterisation and validation of an androgen receptor (AR) transactivation assay for detection of AR agonists and antagonists using a stably transfected human prostate cancer cell line. This 22Rv1/mouse mammary tumour virus glucocorticoid knock-out cell line based AR transactivation assay was validated by criteria in Organisation for Economic Cooperation and Development Guidance Document 34 to determine if the assay performed equally well to the AR EcoScreen Assay included in Test Guideline for AR Transactivation (OECD TG 458). There was no Glucocorticoid Receptor (GR) crosstalk, and no changes in the AR DNA sequence in cells after the successful knock out of GR. Subsequently, the concordance of classifications of the 22 test chemicals was 100% in all laboratories. The AR agonistic and antagonistic inter-laboratory coefficients of variation based on log[10% effect for 10 nM DHT, PC10] and log[inhibitory response of 800 pM DHT by at 30%, IC30] from comprehensive tests were 2.75% and 2.44%, respectively. The AR agonist/antagonist test chemical classifications were consistent across AR EcoScreen ARTA assay data for 82/89%, and the balanced accuracy, sensitivity, and specificity were 83/90%, 88/100% and 78/80%, respectively. This assay was successfully validated and was approved for inclusion in TG 458 in 2020.

Exploration of xenobiotic metabolism within cell lines used for Tox21 chemical screening.

The Tox21 Program has investigated thousands of chemicals with high-throughput screening assays using cell-based assays to link thousands of chemicals to individual molecular targets/pathways. However, these systems have been widely criticized for their suspected lack of 'metabolic competence' to bioactivate or detoxify chemical exposures. In this study, 9 cell line backgrounds used in Tox21 assays (i.e., HepG2, HEK293, Hela, HCT116, ME180, CHO-K1, GH3.TRE-Luc, C3H10T1/2 and MCF7) were evaluated via metabolite formation rates, along with metabolic clearance and metabolite profiling for HepG2, HEK293, and MCF-7aroERE, in comparison to pooled donor (50) suspensions of primary human hepatocytes (PHHs). Using prototype clinical drug substrates for CYP1A2, CYP2B6, and CYP3A4/5, extremely low-to-undetectable CYP450 metabolism was observed (24 h), and consistent with their purported 'lack' of metabolic competence. However, for Phase II metabolizing enzymes and metabolic clearance, surprisingly proficient metabolism was observed for bisphenol AF, bisphenol S, and 7-hydroxycoumarin. Here, comparatively low glucuronidation relative to sulfation was observed in contrast to equivalent levels in PHHs. Overall, while a lack of CYP450 metabolism was confirmed in this benchmarking effort, Tox21 cell lines were not 'incompetent' for xenobiotic metabolism, and displayed surprisingly high proficiency for sulfation that rivaled PHHs. These findings have implications for the interpretation of Tox21 assay data, and establish a framework for evaluating of 'metabolic competence' with in vitro models.

High-throughput toxicogenomic screening of chemicals in the environment using metabolically competent hepatic cell cultures.

The ToxCast in vitro screening program has provided concentration-response bioactivity data across more than a thousand assay endpoints for thousands of chemicals found in our environment and commerce. However, most ToxCast screening assays have evaluated individual biological targets in cancer cell lines lacking integrated physiological functionality (such as receptor signaling, metabolism). We evaluated differentiated HepaRGTM cells, a human liver-derived cell model understood to effectively model physiologically relevant hepatic signaling. Expression of 93 gene transcripts was measured by quantitative polymerase chain reaction using Fluidigm 96.96 dynamic arrays in response to 1060 chemicals tested in eight-point concentration-response. A Bayesian framework quantitatively modeled chemical-induced changes in gene expression via six transcription factors including: aryl hydrocarbon receptor, constitutive androstane receptor, pregnane X receptor, farnesoid X receptor, androgen receptor, and peroxisome proliferator-activated receptor alpha. For these chemicals the network model translates transcriptomic data into Bayesian inferences about molecular targets known to activate toxicological adverse outcome pathways. These data also provide new insights into the molecular signaling network of HepaRGTM cell cultures.

Evaluation of a multiplexed, multispecies nuclear receptor assay for chemical hazard assessment.

Sensitivity to potential endocrine disrupting chemicals in the environment varies across species and is influenced by sequence conservation of their nuclear receptor targets. Here, we evaluated a multiplexed, in vitro assay testing receptors relevant to endocrine and metabolic disruption from five species. The TRANS-FACTORIAL™ system of human nuclear receptors was modified to include additional species: mouse (Mus musculus), frog (Xenopus laevis), zebrafish (Danio rerio), chicken (Gallus gallus), and turtle (Chrysemys picta). Receptors regulating endocrine function and xenobiotic recognition were included, specifically: ERα, ERß, AR, TRα, TRß, PPARγ and PXR. The assay, ECOTOX-FACTORIAL™, was evaluated with 191 chemicals enriched with known receptor ligands. Hierarchical clustering of potency values demonstrated strong coherence of receptor families. Interspecies comparisons of responses within a receptor family showed moderate to high concordance for potencies under 50 µM. PPARγ showed high concordance between mammalian species, 89%, but only 63% between mammalian and zebrafish. For chemicals with potencies below 1 µM, concordances were 89-100% for all receptors except PXR. Concordance showed a strong positive relationship to ligand-binding domain sequence similarity and critical amino acid residues obtained by the Sequence Alignment to Predict Across Species Susceptibility (SeqAPASS) tool. In combination with SeqAPASS, ECOTOX-FACTORIAL may provide efficient screening of important receptors to identify species of high priority for effects monitoring.

The Tox21 10K Compound Library: Collaborative Chemistry Advancing Toxicology.

Since 2009, the Tox21 project has screened ∼8500 chemicals in more than 70 high-throughput assays, generating upward of 100 million data points, with all data publicly available through partner websites at the United States Environmental Protection Agency (EPA), National Center for Advancing Translational Sciences (NCATS), and National Toxicology Program (NTP). Underpinning this public effort is the largest compound library ever constructed specifically for improving understanding of the chemical basis of toxicity across research and regulatory domains. Each Tox21 federal partner brought specialized resources and capabilities to the partnership, including three approximately equal-sized compound libraries. All Tox21 data generated to date have resulted from a confluence of ideas, technologies, and expertise used to design, screen, and analyze the Tox21 10K library. The different programmatic objectives of the partners led to three distinct, overlapping compound libraries that, when combined, not only covered a diversity of chemical structures, use-categories, and properties but also incorporated many types of compound replicates. The history of development of the Tox21 "10K" chemical library and data workflows implemented to ensure quality chemical annotations and allow for various reproducibility assessments are described. Cheminformatics profiling demonstrates how the three partner libraries complement one another to expand the reach of each individual library, as reflected in coverage of regulatory lists, predicted toxicity end points, and physicochemical properties. ToxPrint chemotypes (CTs) and enrichment approaches further demonstrate how the combined partner libraries amplify structure-activity patterns that would otherwise not be detected. Finally, CT enrichments are used to probe global patterns of activity in combined ToxCast and Tox21 activity data sets relative to test-set size and chemical versus biological end point diversity, illustrating the power of CT approaches to discern patterns in chemical-activity data sets. These results support a central premise of the Tox21 program: A collaborative merging of programmatically distinct compound libraries would yield greater rewards than could be achieved separately.

Harmonized Cross-Species Assessment of Endocrine and Metabolic Disruptors by Ecotox FACTORIAL Assay.

Environmental pollution is a threat to humans and wildlife species. Of particular concern are endocrine disrupting chemicals (EDCs). An important target of EDCs is nuclear receptors (NRs) that control endocrine and metabolic responses through transcriptional regulation. Owing in part to structural differences of NRs, adverse effects of EDCs vary significantly among species. Here, we describe a multiplexed reporter assay (the Ecotox FACTORIAL) enabling parallel assessment of compounds' effects on estrogen, androgen, thyroid, and PPARγ receptors of representative mammals, birds, reptiles, amphibians, and fish. The Ecotox FACTORIAL is a single-well assay comprising a set of species-specific, one-hybrid GAL4-NR reporter constructs transiently transfected into test cells. To harmonize cross-species assessments, we used a combination of two approaches. First, we used the same type of test cells for all reporters; second, we implemented a parallel detection of reporter RNAs. The assay demonstrated excellent quality, reproducibility, and insignificant intra-assay variability. Importantly, the EC50 values for NR ligands were consistent with those reported for conventional assays. Using the assay allowed ranking the hazard potential of environmental pollutants (e.g., bisphenols, polycyclic aromatic hydrocarbons, and synthetic progestins) across species. Furthermore, the assay permitted detecting taxa-specific effects of surface water samples. Therefore, the Ecotox FACTORIAL enables harmonized assessment of the endocrine and metabolic disrupting activity of chemicals and surface water in humans as well as in wildlife species.

Selecting a minimal set of androgen receptor assays for screening chemicals.

Screening certain environmental chemicals for their ability to interact with endocrine targets, including the androgen receptor (AR), is an important global concern. We previously developed a model using a battery of eleven in vitro AR assays to predict in vivo AR activity. Here we describe a revised mathematical modeling approach that also incorporates data from newly available assays and demonstrate that subsets of assays can provide close to the same level of predictivity. These subset models are evaluated against the full model using 1820 chemicals, as well as in vitro and in vivo reference chemicals from the literature. Agonist batteries of as few as six assays and antagonist batteries of as few as five assays can yield balanced accuracies of 95% or better relative to the full model. Balanced accuracy for predicting reference chemicals is 100%. An approach is outlined for researchers to develop their own subset batteries to accurately detect AR activity using assays that map to the pathway of key molecular and cellular events involved in chemical-mediated AR activation and transcriptional activity. This work indicates in vitro bioactivity and in silico predictions that map to the AR pathway could be used in an integrated approach to testing and assessment for identifying chemicals that interact directly with the mammalian AR.

Tox21BodyMap: a webtool to map chemical effects on the human body.

To support rapid chemical toxicity assessment and mechanistic hypothesis generation, here we present an intuitive webtool allowing a user to identify target organs in the human body where a substance is estimated to be more likely to produce effects. This tool, called Tox21BodyMap, incorporates results of 9,270 chemicals tested in the United States federal Tox21 research consortium in 971 high-throughput screening (HTS) assays whose targets were mapped onto human organs using organ-specific gene expression data. Via Tox21BodyMap's interactive tools, users can visualize chemical target specificity by organ system, and implement different filtering criteria by changing gene expression thresholds and activity concentration parameters. Dynamic network representations, data tables, and plots with comprehensive activity summaries across all Tox21 HTS assay targets provide an overall picture of chemical bioactivity. Tox21BodyMap webserver is available at https://sandbox.ntp.niehs.nih.gov/bodymap/.

Nontarget Screening of Per- and Polyfluoroalkyl Substances Binding to Human Liver Fatty Acid Binding Protein.

More than 1000 per- and polyfluoroalkyl substances (PFASs) have been discovered by nontarget analysis (NTA), but their prioritization for health concerns is challenging. We developed a method by incorporating size-exclusion column co-elution (SECC) and NTA, to screen PFASs binding to human liver fatty acid binding protein (hL-FABP). Of 74 PFASs assessed, 20 were identified as hL-FABP ligands in which eight of them have high binding affinities. Increased PFAS binding affinities correlate with stronger responses in electrospray ionization (ESI-) and longer retention times on a C18 column. This is well explained by a mechanistic model, which revealed that both polar and hydrophobic interactions are crucial for binding affinities. Encouraged by this, we then developed an SECC method to identify hL-FABP ligands, and all eight high-affinity ligands were selectively captured from 74 PFASs. The method was further applied to an aqueous film-forming foam (AFFF) product in which 31 new hL-FABP ligands were identified. Suspect and nontargeted screening revealed these ligands as analogues of perfluorosulfonic acids and homologues of alkyl ether sulfates (C8- and C10/EOn, C8H17(C2H4O)nSO4-, and C10H21(C2H4O)nSO4-). The SECC method was then applied to AFFF-contaminated surface waters. In addition to perfluorooctanesulfonic acid and perfluorohexanesulfonic acid, eight other AFFF chemicals were discovered as novel ligands, including four C14- and C15/EOn. This study implemented a high-throughput method to prioritize PFASs and revealed the existence of many previously unknown hL-FABP ligands.